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We present a benchmark study for the adsorption of a large 7r-conjugated organic molecule on 
different noble metal surfaces, which is based on x-ray standing wave (XSW) measurements and 
density functional theory calculations with van der Waals (vdW) interactions. The bonding distances 
of diindenoperylene on Cu(lll), Ag(lll), and Au(lll) surfaces (2.51, 3.01, and 3.10 A, respectively) 
determined with the normal incidence XSW technique are compared with calculations. Excellent 
agreement with the experimental data, i.e., deviations less than 0.1 A, is achieved using the Perdew- 
Burke-Ernzerhof (PBE) functional with vdW interactions that include the collective response of 
substrate electrons (the PBE+vdW sur method). It is noteworthy that the calculations show that 
the vdW contribution to the adsorption energy increases in the order Au(lll) < Ag(lll) < Cu(lll). 

PACS numbers: 68.49.Uv, 68.43.-h, 71.15.Mb, 87.15.A- 
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I. INTRODUCTION 

The reliable prediction of the equilibrium structure 
and energetics of hybrid inorganic/organic systems from 
first principles represents a great challenge for theoret- 
ical methods due to the interplay of covalent interac- 
tions, electron transfer processes, Pauli repulsion, and 
van der Waals (vdW) interactions. During recent years, 
huge efforts have been made to incorporate vdW interac- 
tions into density functional theory (DFT) calculations 
in order to determine the structure and stability of n- 
conjugated organic molecules on solid surfaces 1 - - — . Un- 
derstanding these interface properties is relevant, inter 
alia, for electron transfer processes in organic devices. 
Until now and despite the obvious benefit, there are 
only few a studies of metal-organic interfaces combin- 
ing theory and experiment. Here, x-ray standing wave 
(XSW) measurements can provide an important test for 
DFT calculations 2 ^. This is particularly important for 
systems with strong vdW contributions to the overall 
bonding, for which no simple substrate dependence is 
expected. 

As model system wc chose diindenoperylene (DIP, 
C32H16), a 7r-conjugatcd organic semiconductor with 
excellent optoelectronic device performance, which has 
been studied over the last decade both in thin-films^— 
and in monolayers on noble metal surfaces^—. With 
respect to its chemical structure, DIP is a relatively 
simple, planar hydrocarbon without heteroatoms. In 
contrast to the intensely studied perylcne derivative 



3,4,9, 10-perylene tetracarboxylic dianhydride (PTCDA, 
C 32 H 8 6 )^^ with its four keto groups, the DIP- 
substrate interaction is not complicated by polar side 
groups, and the influence of intermolecular interactions 
is expected to be smaller than for PTCDA 2 ^. Here, we 
present a systematic study with high-precision experi- 
mental data and state-of-the-art calculations of DIP ad- 
sorbed on Cu(lll), Ag(lll), and Au(lll). This allows 
us to assess the role and relative contribution of the vdW 
interactions, which, contrary to simplistic pictures, we 
find here to be lowest for the most polarizable substrate. 



II. COMPUTATIONAL APPROACH 

DFT calculations were performed using a method that 
extends standard pairwisc vdW approache o 22 i 23 to model 
adsorbatcs on surfaces^. This was achieved by combin- 
ing the DFT+vdW scheme 2 ^ with the Lifshitz-Zaremba- 
Kohn (LZK) theory for vdW interaction between an 
atom and the surfac e) 24 ' 25 of a solid. In our approach 
(DFT+vdW s ) , the vdW energy is given by a sum 
of CQ b R~® terms, where R a b are the distances between 
atoms a and b, in analogy to standard pairwise disper- 
sion corrected DFT methods. However, by employing 
the LZK theory we include the many-body collective re- 
sponse (screening) of the substrate electrons in the de- 
termination of the Cq coefficients and vdW radii, going 
effectively beyond the pairwise description. Interface po- 
larization effects arc accounted for via the inclusion of 



semi-local hybridization due to the dependence of the 
Cg 6 interatomic coefficients on the electron density in 
the DFT+vdW method. The DFT+vdW surf method has 
been shown to yield remarkably accurate results for the 
structure and adsorption energies of xenon, benzene, and 
PTCDA on a variety of (111) metal surfaces^. The 
FHI-aims code2£ was employed for our DFT calculations. 
The repeated-slab method was used to model all systems 
with the vacuum gap set to 20 A. In all calculations, con- 
vergence criteria of 10 -5 electrons for the electron den- 
sity and 10 -6 eV for the total energy of the system were 
used. A convergence criterion of 0.01 eV/A for the max- 
imum final force was used for all structure relaxations. 
The scaled zeroth-order regular approximation (ZORA) 
was applied for inclusion of scalar relativistic effects^. 
The DFT+vdW s method employed the Pcrdcw-Burkc- 
Ernzerhof (PBE) functional^. The sampling of the Bril- 
louin zone was done using a (2 x 2 x 1) fc-point grid. 

We used a (7 x 7) unit cell composed of a metal sur- 
face of three layers and one single DIP molecule. In the 
absence of experimental data for the in-plane registry we 
placed the central ring of the molecule aligned with a top- 
most metal layer atom and the major axis of the molecule 
aligned along the cell diagonal. This structure was 
adopted for the Cu(lll), Ag(lll), and Au(lll) surfaces. 
In each simulation we obtained the adsorption energy 
curve using a rigid DIP molecule and tuning the surface- 
molecule distance d. The adsorption energy per molecule 
£ a d s was calculated from £ ads = E tot - (E sulf + £bip), 
where -E sur f is the energy per unit cell of the isolated 
metal surface, -Edip is the energy per unit cell of the iso- 
lated DIP molecule, and E to t is the energy per unit cell 
of the combined system. We also obtained the relaxed 
geometries for all three systems starting from the static 
equilibrium geometry. During geometry relaxation, we 
allowed only the topmost metal layer and the molecule 
to relax while the other two metal layers were fixed. From 
the final relaxed configurations we obtained the bonding 
distance d by taking the average position of all DIP car- 
bon atoms with respect to the unrelaxed topmost surface 
layer. This definition is consistent with the analysis of the 
XSW data. 



III. EXPERIMENTAL DETAILS 

To measure the bonding distance of DIP we used 
the XSW technique^, which yields precise and element- 
specific structural data. The experiments were performed 
at beamlinc ID32 of the ESRF^i. DIP films were pre- 
pared and studied in situ under ultrahigh-vacuum con- 
ditions. A separate preparation chamber contained a 
Knudsen cell, a quartz crystal microbalance, installations 
for Ar + sputtering, and a temperature-controlled sam- 
ple stage. The main chamber, in which the XSW mea- 
surements were performed, was equipped with a sample 
manipulator and a hemispherical SPECS PHOIBOS 225 
HV photoclcctron analyzer. The XSW experiments were 
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FIG. 1: (Color online) C Is core-level shift observed for a sub- 
monolayer of DIP (inset) on Cu(lll), Ag(lll), and Au(lll). 
The spectra were taken at an emission angle of 45° with 
the XSW setup at ID32. From each signal a Shirley back- 
ground was subtracted and then fitted with a Voigt function 
for the main peak and a Gaussian function for possible shake- 
up peaks. 



carried out at room temperature in back-reflection geom- 
etry using the (111) Bragg reflection of the crystals for at 
least two films per substrate to check for reproducibility 
of the results. The detection angle of the analyzer was 
~90° relative to the surface normal with an acceptance 
angle of ±7.5°. We note that in this configuration non- 
dipolar contributions to the photoelecton yield can be ef- 
fectively avoided^. The Cu(lll), Ag(lll), and Au(lll) 
single crystals were mounted on different sample hold- 
ers for individual treatment. The surfaces were prepared 
by repeated cycles of Ar + bombardment and annealing 
at 700 K. Surface cleanliness was confirmed with x-ray 
photoelectron spectroscopy (XPS) as well as low-energy 
electron diffraction (LEED). Sublimation grade DIP was 
evaporated from a home-built Knudsen cell. The inten- 
sity ratio of the C Is signal relative to a substrate core- 
level, normalized with the corresponding photocmission 
cross sections, was used to determine the number of DIP 
molecules on the surface. With the unit cell size of DIP 
on Cu(lll)i2, Ag(lll)^, and Au(lll)ii, the coverages 
were calculated to be between 0.3 and 0.9 ML. 



IV. RESULTS AND ANALYSIS 

A. Experimental results 

The C Is core-level signals of DIP on Cu(lll), 
Ag(lll), and Au(lll), which were used for the XSW 
measurements, are shown in Fig. [1] The main peaks are 
expected to consist of two principal components (C-C vs 
C-H bound atoms) which, however, could not be resolved 
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FIG. 2: (Color online) (a) Typical XSW data for DIP show- 
ing the reflectivity (triangles) and photoelectron yield (cir- 
cles) of the C Is signal on Cu(lll), Ag(lll), and Au(lll). 
The solid lines correspond to least-squares fits of the re- 
flectivity and photoelectron yield, which reveal the coherent 
position Ph and coherent fraction fn- Bragg energies are 
£ B ragg = 2.97 keV [Cu(lll)] and £ B rag g = 2.63 keV [Ag(lll) 
and Au(lll)]. (b) All XSW results for DIP on Cu(lll), 
Ag(lll), and Au(lll) displayed in an Argand diagram. Here, 
each datapoint represents one single XSW measurement yield- 
ing fg (length of a vector) and Ph (angle of a vector). The 
three vectors point to the average values of fn and Ph for 
DIP on Cu(lll), Ag(lll), and Au(lll). Film 1 of DIP on 
Ag(lll) was measured with a different crystal compared to 
films 2 and 3. 



with the energy resolution of the XSW setup. In addi- 
tion to each main peak, a second weak feature at ~1 cV 
higher binding energy possibly related to a shake-up pro- 
cess can be observed. Obviously, the binding energy of 
the C Is main line of DIP follows E^ g > Eg u > £ Au , be- 
ing 284.5 eV on Ag(lll), 284.2 eV on Cu(lll), and 283.7 
cV on Au(lll). Furthermore, the C Is peak of DIP on 
Ag(lll) exhibits a stronger asymmetry than on Cu(lll) 
and Au(lll)i^ A detailed discussion of the spectroscopic 



TABLE I: Results of XSW experiments: Coherent fraction 
fn, coherent position Ph, and bonding distance dn of DIP on 
the three noble metals. The parameters refer to an average of 
several XSW measurements with the corresponding standard 
deviation as error bars. 





Ih 


Ph 


d H 


Cu(lll) 
Ag(lll) 
Au(lll) 


0.48 ± 0.09 
0.55 ± 0.08 
0.62 ±0.06 


0.20 ±0.01 
0.28 ±0.02 
0.35 ±0.01 


(2.51 ±0.03) A 
(3.01 ±0.04) A 
(3.17 ±0.03) A" 



"By taking the surface reconstruction of Au(lll) into account, 
dn is reduced to 3.10 A. 



features is beyond the scope of this paper in which we fo- 
cus on the XSW results. 

Representative results of the XSW experiments are 
shown in Fig. [2ja). In each panel the measured reflec- 
tivity of the substrate and the corresponding C 1 s pho- 
toelectron yield is displayed. Least-squares fits of the 
data give the coherent position Ph and hence the aver- 
age bonding distance dn — do(l + Ph)^, where do is the 
substrate lattice plane spacing. Based on results of all 
XSW experiments we calculate the average bonding dis- 
tance dn and the standard deviation; see Fig. [U[b). For 
Cu(lll) we thus find (2.51 ± 0.03) A, and for Ag(lll) 
(3.01 ±0.04) A. Due to the reconstruction of the Au(lll) 
surface, which results in a 3% larger spacing between the 
first and second Au layers^, the bonding distance de- 
creases from the measured apparent value (3.17 ± 0.03) 
A to (3.10 ± 0.03) A. All experimental results are sum- 
marized in Table HI Note that although the coverage of 
the two (three) DIP films prepared on each substrate was 
not identical, we did not observe a significantly coverage- 
dependent bonding distance dn- 

Comparing these results with the bonding distances 
of PTCDA on the same metal surfaces, i.e., dn = 2.66 
A on Cu(lll)A d H = 2.86 A on Ag(lll)^i£, and 
dn = 3.27 A on Au(lll)jil we see that the bonding 
distances follow the same order, i.e., dn (Cu) < dfj(Ag) 
< dy(Au). Moreover, the results demonstrate that the 
absence of the C=0 groups in DIP affects the bonding 
distance of the molecule only weakly. 



B. Computational results 

Having established precise experimental data, we now 
turn to the results of our DFT calculations. The aver- 
age bonding distances of DIP obtained from fully relaxed 
structures arc d = 2.59A on Cu(lll), d = 2.94A on 
Ag(lll), and d = 3.22 A on Au(lll), sec Table M and 
Fig. [3] We hence find that the PBE+vdW surf method 
applied to DIP on Cu(lll), Ag(lll), and Au(lll) yields 
an agreement better than 0.1 A between theoretical cal- 
culations and experiments. In accordance with the bond- 
ing distances, the calculated adsorption energies listed 
in Table |TT] follow the trend |£: ads (Cu)| > |£; ads (Ag)| > 



TABLE II: Adsorption energy -E a d s of the relaxed structures, 
vdW 8urf binding energy -E v dw in parentheses as derived from 
data shown in Fig. [3] and distances d between the topmost 
layer of the metal and the carbon backbone of DIP. d m i n / max 
refer to the lowest/highest bonding distances of a carbon atom 
within a DIP molecule. 



DISCUSSION 





-Bads (-B v dw) 


d 


Gfmin 


0/xaax 


Cu(lll) 
Ag(lll) 
Au(lll) 


-4.74 (-5.28) eV 
-3.55 (-4.56) eV 
-2.53 (-3.06) eV 


2.59 A 
2.94 A 
3.22 A 


2.38 A 
2.89 A 
3.15 A 


2.79 A 
3.01 A 
3.29 A 



|-Sads(Au)|. Interestingly, Fig. [3] shows that on Cu(lll) 
the Pauli repulsion sets in rather weakly [a less steep 
E a ds(d) for small distances] compared to Ag(lll) and 
Au(lll), which is due to significant interaction between 
DIP and Cu(lll). One may speculate that the interac- 
tion mechanism includes hybridization between DIP and 
Cu states. 

In addition to the adsorption energies and average 
bonding distances, Table [TT1 holds the minimal and max- 
imal bonding distances rf m i n / max of individual carbon 
atoms in DIP. These values indicate that the molecule 
adsorbs in a slightly tilted or distorted geometry. For 
Cu(lll), where the effect is most pronounced, the calcu- 
lated bonding distances d m i n and d max differ by ~0.4A, 
which is equivalent to a molecular tilt angle of 1.5°. The 
corresponding spread of vertical positions of the carbon 
atoms leads to a reduced fjj in the XSW scans. Model 
simulations similar to those presented in Ref. [3J show 
that the DFT-derived adsorption geometry on Cu(lll) 
results in a relatively small decrease of the coherent frac- 
tion (Afjj = —0.07), which lies within the standard de- 
viation of our XSW measurements. 

To obtain a better understanding of the influence of 
lateral intermolecular interactions on the DIP adsorp- 
tion geometry, we also computed the relaxed DIP geom- 
etry for different Cu(lll), Ag(lll), and Au(lll) unit 
cells. For DIP on Cu(lll), we increased the unit cell 
from (7 x 7) to (9 x 7) in order to reduce the molecule- 
molecule interactions. We studied various configurations, 
finding a flat relaxed geometry for each case considered. 
The bonding distance is slightly larger (2.64 A) than for 
the calculation with the smaller unit cell. For DIP on 
Ag(lll), we also considered a unit cell which was de- 
termined from a closed packed monolayer on Ag(lll)^. 
The relaxed geometry of the molecules in the monolayer 
is flat and the bonding distance d = 2.99 A in almost 
perfect agreement with the experimental one, i.e., even 
better than the result for Ag(lll) shown in Table ITT1 For 
a (9 x 5) unit cell of Au(lll), the relaxed DIP geometry 
yields an equilibrium distance of 3.15 A, also in slightly 
better agreement with experiment than the result shown 
in Table HU Overall, these calculations agree with the 
experimental observation that the vertical DIP position 
depends only weakly on surface coverage. 



O DIP/Cu(lll)-PBE 
A DIP/Ag(lll)-PBE 
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FIG. 3: (Color online) Adsorption energy 25 a ds for the un- 
relaxed DIP molecule as a function of its averaged distance 
d from the Cu(lll), Ag(lll), and Au(lll) surfaces. The 
curves are shown for the PBE functional with and without the 
inclusion of long-range vdW interactions using the vdW sur 
method. The reported contribution of the vdW energy is 
shown at the equilibrium distance corresponding to the fully 
relaxed DIP-surface geometry (see text). 

With the experimental and theoretical values at hand, 
and in view of their excellent agreement, we are in a good 
position to discuss the vdW interactions and the bonding 
distances in more detail. As described above, the (atom- 
atom) vdW energy is computed as Cg b R~^, where the 
Cg fe coefficient determines the strength of the interac- 
tion between atoms a and b, while R a b is the distance 
between adsorbate and substrate atoms (Fig. 0]). Inte- 
gration of the vdW energy for a single atom adsorbed on 
a semi-infinite surface yields the atom-surface vdW en- 
ergy a o 35 ' 36 C^~ s (z — zo)~ 3 , where now C^' s determines 
the interaction strength between atom and surface, z cor- 
responds to the distance of the atom to the uppermost 
surface layer, and zq indicates the position of the surface 
image plane. In a rather naive picture, the C^ coef- 
ficient can be determined simply from the Cg a and the 



Cg coefficients that correspond to the adsorbed atom 
and the metal atom, respectively. However, the situation 
for real surfaces is more complex because both localized 
and bulk metal electrons contribute to the C^~ s coeffi- 
cient in a non-trivial way, meaning that this coefficient 
depends on the dielectric function of the underlying solid. 
We computed the C^~ s coefficients corresponding to the 
interaction between a carbon atom and the Cu(lll), 
Ag(lll), and Au(lll) surfaces. When describing the 
metal surface as a simple collection of non-interacting 
atoms we obtain Cl?" Cu 



0.68, Cf~ Ag = 0.55, and 



C^' Au = 0.50 hartree-bohr 3 . In contrast, when using the 
more appropriate Lifshitz-Zarcmba-Kohn cxpressio n 24 i 25 



for C^~ s , we 

£iC-Au 



obtain C^ Cu 



0.35, a 



C- 



0.35, and 



0.33 hartree-bohr 3 . This clearly illustrates that 
the vdW interaction between an atom and a solid sur- 
face is significantly modified by the collective electronic 
response within the substrate surface£ ' 35 i 36 . The very 
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ever, at shorter molecule-surface distances, which include 
the equilibrium distance, the adsorption energy is de- 
termined by an interplay between the vdW attraction 
and the Pauli repulsion with a possible covalent compo- 
nent. The Pauli repulsion follows roughly the trend of 
decreasing vdW radii, with a faster onset in terms of the 
molecule-surface distance for Au (with the largest vdW 
radius), and then decreasing for Ag and Cu. Therefore, 
for Au the balance between the vdW attraction and the 
Pauli repulsion is obtained further away from the sub- 
strate (i.e., at larger adsorption distances) than for Cu, 
which in turn makes the adsorption energies lower for Au 
than for Cu, in contrast to the possible naive expectation 
of Au with its higher polarizability and Cq coefficient ex- 
hibiting a stronger vdW interaction than Cu. 

The difference in the vdW energy distribution for 
DIP on Cu(lll), Ag(lll), and Au(lll) is visualized in 
Fig. |U[b), where the vdW energy between DIP and sub- 
strate atoms is plotted as a function of their distance Rab- 
in contrast to Ag(lll) and Au(lll), the small bonding 
distance of DIP on Cu(lll) results in a second peak in the 
histogram at ~3.6 A, which originates from the higher 
atomic density of the Cu substrate. 



VI. CONCLUSION 

In conclusion, the bonding distances calculated with 
the PBE+vdW surf method are in excellent agreement 
with the XSW data for DIP on Cu(lll), Ag(lll), and 
Au(lll) (2.51, 3.01, and 3.10 A, respectively). Our 
combined study demonstrates that the vdW energy is 
larger for DIP on Cu(lll) than for DIP on Ag(lll) and 
Au(lll). Future investigations on the electronic prop- 
erties of these systems, which can draw on the findings 
presented here, will contribute to an even better under- 
standing of the adsorption process. 



FIG. 4: (Color online) (a) Schematic picture of the adsorp- 
tion of DIP on a (111) crystal. Each arrow corresponds to a 
certain distance R a t, between adsorbate and substrate atoms, 
(b) vdW energy of DIP adsorbed on Cu(lll), Ag(lll), and 
Au(lll) as a function of the distance between adsorbate and 
substrate atoms. 

similar LZK C3 coefficients for Cu, Ag, and Au lead to es- 
sentially the same adsorption energy at large distances for 
DIP on Cu(lll), Ag(lll), and Au(lll) (Fig.©. How- 
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